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ABSTRACT 

(-H \ Ultraviolet and X-ray observations show evidence of outflowing gas around many active 

pi J galactic nuclei. It has been proposed that some of these outflows are driven off gas infalling 

towards the central supermassive black hole. We perform radiative transfer calculations to 
compute the gas ionization state and the emergent X-ray spectra for both two- and three- 
dimensional (3D) hydrodynamical simulations of this outflow-from-inflow scenario. By com- 
parison with observations, our results can be used to test the theoretical models and guide fu- 
ture numerical simulations. We predict both absorption and emission features, most of which 
are formed in a polar funnel of relatively dense (10^^" - 10~^^ g cm~^) outflowing gas. This 
outflow causes strong absorption for observer orientation angles of ^ 35 degrees. Particularly 
in 3D, the strength of this absorption varies significantly for different lines-of-sight owing to 
the fragmentary structure of the gas flow. Although infalling material occupies a large fraction 
of the simulation volume, we do not find that it imprints strong absorption features in the X- 
ray spectra since the ionization state is predicted to be very high. Thus, an absence of observed 
inflow absorption features does not exclude the models. The main spectroscopic consequence 
of the infalling gas is a Compton-scattered continuum component that partially re-fills the 
absorption features caused by the outflowing polar funnel. Fluorescence and scattering in the 
outflow is predicted to give rise to several emission features including a multi-component 
Fe Ka emission complex for all observer orientations. For the hydrodynamical simulations 
considered, we predict both ionization states and column densities for the outflowing gas 
that are too high to be quantitatively consistent with well-observed X-ray absorption systems. 
Nevertheless, our results are qualitatively encouraging and further exploration of the model 
parameter space is warranted. Higher resolution hydrodynamic simulations are needed to de- 
termine whether the outflows fragment on scales unresolved in our current study, which may 
yield the denser lower ionization material that could reconcile the models and the observa- 
tions. 

Key words: hydrodynamics - radiative transfer - methods: numerical - galaxies: active - 
X-rays: galaxies 



1 INTRODUCTION 

The energy liberated during accretion by supermassive black holes 
can be sufficient that active galactic nuclei (AGN) could signifi- 
cantly influence their environments. Indeed, it has been argued that 
feedback from AGN may play an important role in the evolution 
of galaxies (e.g. Croton et al. 2006). One potential mechanism by 
which AGN may affect their host galaxy is via energy and mo- 
mentum carried by winds or outflows (King 2003, 2005; Ciotti & 
Ostriker 2007; Ciotti et al. 2009, 2010; Ostriker et al. 2010). 



There is copious observational evidence suggesting that many 
AGN have associated outflows. In particular, blue-shifted absorp- 
tion lines are commonly observed in both ultraviolet (uv) and X- 
ray spectra (see reviews by Crenshaw et al. 2003; Turner & Miller 
2009). These absorption lines suggest outflows exist with a wide 
variety of properties: although the shifts of many uv and soft X-ray 
lines typically suggest outflow speeds of only a few hundreds to 
thousands of km s~^, uv spectra of broad absorption line quasars 
show line shifts up to ~ 0.1c (Weymann et al. 1981) and evidence 
of even higher velocity flows has been found in X-ray observations 
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(e.g. Reeves et al. 2003; Chartas et al. 2002; Tombesi et al. 2010). 
In addition, detection of rapid (~ ks) variability in some highly ion- 
ized absorption components implies that that they are located rela- 
tively close (~ hundreds of gravitational radii) to the AGN (Reeves 
et al. 2004; Miller et al. 2007; Turner et al. 2007, 2011). In other 
cases, however, lack of variability has been used to argue that com- 
ponents of the X-ray absorbing gas are distant from the AGN, likely 
located on parsec scales from the central black hole (Netzer et al. 
2003; Behar et al. 2003; Krongold et al. 2010; Kaastra et al. 2012). 
This wide diversity of absorber properties suggests that there may 
be multiple means by which AGN drive outflows and so theoretical 
modelling for a variety of outflow mechanisms is required to build 
a complete picture. 

Since AGN are extremely luminous, it is natural to investigate 
the role of radiation in ejecting material from their environment. In 
a previous study (Sim et al. 2010b), we computed synthetic X-ray 
spectra for simulations of a line-driven AGN accretion disk wind 
(Proga & Kallman 2004) launched from within a few hundred grav- 
itational radii of a supermassive black hole. We found that such 
an outflow might account for a variety of observable spectral fea- 
tures associated with highly ionized, high- velocity gas. However, 
many observed spectral features (e.g. absorption lines with rela- 
tively small blueshifts or narrow components of emission lines) are 
unlikely to be formed in outflows very close to the central black 
hole. These are more probably associated with distant material 
where characteristic flow velocities are expected to be smaller. One 
promising origin for these distant outflows is a wind blown from 
the inner regions of a torus of dusty material surrounding the AGN 
(Balsara & Krolik 1993; Dorodnitsyn et al. 2008a,b; Dorodnitsyn 
& Kallman 2009). In that model, the mass reservoir for the outflow 
is rotationally supported cold material and the resulting outflow has 
velocities of ~ 1000 km s~^ close to the rotation axis, dropping to 
a few hundred km s^^ at higher inclinations (Dorodnitsyn et al. 
2008b). Such a model predicts observable outflow features which 
are broadly compatible with the properties of many observed X-ray 
absorbers. 

An alternative model is that in which an outflow is driven not 
from a rotationally supported torus but directly from gas infalling 
towards a supermassive black hole. This scenario has been mod- 
elled on parsec scales by Proga (2007), Proga et al. (2008), Kuro- 
sawa & Proga (2009b, hereafter KP09) and Kurosawa & Proga 
(2009a). It was found that momentum and energy supplied by radi- 
ation from the AGN could strongly affect the dynamics of infalling 
gas, giving rise to a characteristic flow pattern which includes a 
polar outflow driven off the infalling gas. Proga (2007) and Proga 
et al. (2008) studied this process using two-dimensional (2D) sim- 
ulations both with and without gas rotation and stressed that the 
formation of such a flow pattern was probable if the AGN luminos- 
ity is larger than around one per cent of the Eddington luminosity. 
KP09 extended the studies to included fully three-dimensional (3D) 
simulations and confirmed that the characteristic flow pattern of the 
earlier 2D calculations was also robust in 3D. Kurosawa & Proga 
(2009a) examined the consequences of coupling the black hole ac- 
cretion luminosity to the mass inflow rate in the simulations and 
still found the same flow pattern for a outflow radiatively driven off 
hot infalling gas in many cases. 

The goal of this paper is to compute synthetic X-ray spectra by 
carrying out radiative transfer simulations based on the KP09 hy- 
drodynamical calculations. These synthetic spectra will allow us to 
quantify the observable consequences of the outflow/inflow pattern 
and asses how well they might account for observed spectroscopic 
features. One key motivation of this work is the question of whether 



the infalling gas in the simulation strongly affects the X-ray spec- 
tra. Signatures of infall are not commonly observed and so it is 
important to consider whether this is compatible with models for 
outflows that are driven off a mass reservoir of infalling gas. We 
also wish to compare the observational features predicted by 3D 
outflow simulations with those of axi-symmetric 2D calculations. 
KP09 showed that the overall flow pattern in 3D was similar to that 
in 2D but with significantly more structure in the polar outflow. Our 
radiative transfer simulations allow us to quantity the observable 
consequences of this extra complexity. 

We begin by summarising the key properties of the KP09 sim- 
ulations and our radiative transfer calculations in Section 2. We 
then present our synthetic spectra in Section 3 and discuss their in- 
terpretation in Section 4. Our conclusions are summarised in Sec- 
tion 5. 



2 NUMERICAL SIMULATIONS 

2.1 Input from radiation hydrodynamical simulations 

As input for our radiative transfer calculations, we have taken 
snapshots of density and velocity from two of the radiation hy- 
drodynamical simulations presented by KP09. Their simulations 
are for the region 3000 < r/r^ < 1.5 x 10*" (where r-g = 
GMbh/c^ = 1.5 X 10^^ cm) around a central black hole of 
mass A/bh ~ 10* Af©. An accretion disk surrounding the black 
hole is assumed to radiate at a bolometric luminosity of Lboi ~ 
7.6 X 10*^ ergs s^^. The simulations were initialized with a uni- 
form mass density (po ~ 10^^^ g cm^^) and gas temperature 
(To = 2 X 10^ K). Radial velocities were initialized to zero but ro- 
tational motions were included in the setup as described in section 
2.3 of KP09. The polar axis is treated with an axis-of-symmetry 
boundary condition while outflow boundary conditions are applied 
at both the inner and outer radial boundaries. The density and tem- 
perature of material entering through the outer boundary of the sim- 
ulation (at r = 1.5 X lO^'r-g ~ 7 pc) were fixed to po and To while 
the radial velocity at the boundary was allowed to float (see KP09 
for full details of the simulation setup^). 

We have studied both the KP09 simulations that included gas 
rotation: their axi-symmetric simulations (model II) and full three- 
dimensional simulation (model IV). The initial conditions for these 
two hydrodynamical simulations were identical - the only dif- 
ference was the assumed axi-symmetry in model II. Thus, com- 
parison of these simulations will allow us to quantify how the 
imposed axi-symmetry of two-dimensional models influences our 
conclusions. We will refer to these models as KP09-II and KP09- 
IV, respectively. In physical units, the time between the snap- 
shots recorded during the KP09 hydrodynamical simulations was 
~ 300 yr. This is too long for it to be useful to compare multi- 
ple snapshots to estimate how the spectral features may vary over 
observable timescales. Consequently we consider only one repre- 
sentative snapshot for each simulation. The particular snapshots we 
adopt are taken from the phase where the hydrodynamical simula- 
tion has evolved into a time-averaged steady state and are the same 
as those shown in figure 2 of KP09. Neither the hydrodynamical 
simulations nor our radiative transfer calculations assume symme- 
try about the xy-plane (i.e. the full interval for the spherical polar 
angle is included, < 9 < 180 deg.). However, the flow structure 

^ For simulation movies, see 

http://www.astro.comell.edu/~kurosawa/research/agn.html. 
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Figure 1. Density distributions in the +z region for tlie models we consider. Tlie top left panel show a slice in the x = 0-plane for our two-dimensional model 
(KP09-II), which is symmetric under rotation about the z-axis. The remaining three panels show slices through the density distribution of our three-dimensional 
model (KP09-IV) in the planes defined by the spherical polar coordinate 4> = 28, 62 and 90 degrees. 



is very similar in the +z and —2 hemispheres and so we will focus 
our discussion on only the +z hemisphere (0 < ^ < 90 deg); in all 
important respects, the —z region mirrors the same behaviour. 

Figure 1 shows two-dimensional slices of the density in the 
model snapshots we have studied (for KP09-IV we include three 
representative slices). In general, these all show similar features, 
as described by KP09. The flow structure can roughly be divided 
into three regimes. Very close to the rotation axis {6 < 15 deg), 
there is a region of very low density gas (p < 10"^^ g cm^'^), 
which is flowing away from the central source. Surrounding this 
is an outflowing funnel of relatively dense gas (typically p > 
2 X 10"^" g cm~^) which is being accelerated out of the system 
(the total mass outflow rate through the outer boundary of the sim- 
ulation is ~ 0.8 M0 yr^^). The gas in the polar funnel is ther- 
mally unstable. The hot thermally stable infalling gas in the funnel 



is brought to the regime of thermal instability by the action of the 
centrifugal force and radiation pressure on electrons - both forces 
slow down the gas and lead to a density increase and fast cool- 
ing. This in turn enhances the outward line force that is sensitive 
to the density, photoionization parameter and optical depth. Con- 
sequently, a fragmented outflow develops with dense substructures 
(exceeding p — 10^^** g cm^'* in places). Outside this outflowing 
stream (9 > 35 deg), the gas density is significantly lower (in most 
of this region, p < lO^'' g cm~"^) and the predominant motion is 
infall towards the central source (the net inflow rate through the 
inner boundary is ~ 0.8 Mq yr~^). The infalling gas outside the 
funnel is supersonic, heated mostly by adiabatic compression and 
its temperature is high so that this gas is always thermally stable. 

The outflowing material is sufficiently dense that it can be ex- 
pected to imprint strong signatures on the spectrum of the source. 
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Figure 2. Hydrogen column density measured through the models from the 
origin to an observer at inclination 9 to the z-axis for our models. For the 
three-dimensional KP09-1V model, we plot the column-density averaged 
over all azimuthal angles. The dashed vertical line marks the divide between 
orientations for which the KP09-II column-density weighted mean radial 
velocity is positive (i.e. the mean motion is outflow) and negative (i.e. mean 
inflow). We show only the +z hemisphere - the same behaviour is mirrored 
in the —z hemisphere. 



Figure 3. Hydrogen column density through model KP09-rV as a function 
of azimuthal angle (<f)) for fixed polar angle (6). Values are shown for five 
0-values: 5, 15, 22, 44 and 70 degrees (see labels in the Figure). 



Figure 2 shows the hydrogen column density from the origin 
through our models as a function of inclination-angle (9, mea- 
sured from the 2-axis). For the KP09-IV model we show the col- 
urrm density averaged over all azimuthal angles. In model KP09- 
II, the column density is relatively modest (A^^h < 10^^ cm~^) 
very close to the polar axis but is significant (~ 10^'^ cm~^) for 
most lines-of-sight through the flow. The highest column densities 
( ^ 3 X 10'^'^ cm^^) are for lines-of-sight that pass through the out- 
flowing polar funnel (15 < ^ < 35 deg). We also note that the 
variation of the column-density with orientation is most complex 
in this region. This is caused by the fragmentation of the outflow- 
ing gas - the very highest column densities are for lines-of-sight 
that pass through dense clouds. It is therefore to be expected that 
the most complex spectral features will be imprinted when the cen- 
tral X-ray source is viewed through this gas. 

In model KP09-IV, the density distribution is similar to that 
in KP09-II except that the outflowing funnel is now able to frag- 
ment into three-dimensional structures (i.e. there is no imposed axi- 
symmetry). Consequently, the density becomes a strong function of 
the azimuthal angle in the vicinity of the polar outflow. This is illus- 
trated in Figure 3, which shows the variation of the column-density 
through our three-dimensional model (KP09-IV) as a function of 
azimuthal angle ((j)) for several fixed values of the polar angle (9). 
For moderate and high inclinations {6 > 40 deg), the azimuthal 
variation of the column density is small. At these inclinations, the 
line-of-sight is mostly passing through the large volume of low- 
density, infalling material, which is not strongly fragmented. How- 
ever, for smaller values of 6, the variation of A^h with (/> becomes 
increasingly significant - for S ~ 5 deg., there is almost an order 
of magnitude of variation. 



2.2 Radiative transfer simulations 

We have carried our radiative transfer calculations for the snap- 
shots from the KP09 simulations described above using our multi- 
dimensional Monte Carlo radiative transfer code (Sim 2005; Sim 
et al. 2008, 2010a). Except where noted below, these simulations 
were performed in an identical manner to that described in Sim 
et al. (2010b): we first carry out an iterative sequence of calcula- 
tions to compute the ionization state of the gas assuming a centrally 
concentrated primary X-ray source. We then use our Monte Carlo 
quanta to construct detailed synthetic X-ray spectra for a range of 
observer orientations. 

Table 1 lists the important parameters adopted for the radiative 
transfer calculations (the same radiation parameters were adopted 
for our simulations of both models KP09-II and KP09-IV). Al- 
though our Monte Carlo method does not rely on a pre-defined pho- 
ton frequency-grid, the effective resolution of the synthetic spec- 
trum is limited by Monte Carlo noise to roughly c AE/E — 
1000 km s"'^ in the simulations presented here. 

As in Sim et al. (2010b), we have assumed that the primary X- 
ray source is centrally located and isotropically emits a power-law 
spectrum with photon index F — 2.1. The source X-ray luminosity 
{Lx = 1-5 X 10''* ergs s'^) is chosen so that we match the flux 
assumed by KP09 in the 1-20 keV region. 

Since the circularization radius for the infalling gas chosen by 
KP09 (rc = ISOOrg) is smaller than the radius of the inner bound- 
ary of the simulation (ri = 3000rg), an optically thick disk does 
not form in the computational domain. Consequently, in contrast to 
Sim et al. (2010b), we do not halt the propagation of our Monte 
Carlo quanta if they reach the 2 = plane but allow them to cross 
unhindered. 
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Table 1. Paiameters for the simulations. 
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Parameter 


Adopted value 


Mass of central black hole, M^i^ 




Bolometiic luminosity, Lbol 


7.6 X 10*5 ergs s"! (~ 0.6LEdd) 


X-ray source luminosity (2-10 keV), Lx 


1.5 X 10'*"' ergs s-^ 


X-ray source power-law photon index, F 


2.1 


Range of source photons in simulation 


0.1-511 keV 


Size of primary emission region, Tcr 


8.8 X 10" cm (6rg) 


Inner radius of simulated region, rmin 


4.4 X lO^*^ cm 


Outer radius of simulated region, rmax 


2.2 X 10" cm 


Size of Cartesian Rad. Trans, simulation grid 


2.5 X 10" cm 


3D Cartesian Rad. Trans, grid cells 


256 X 256 X 256 
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Figure 4. Mean charge of Fe ions (upper left) and the kinetic temperature (upper right) obtained from our radiative transfer simulations for the KP09-II model. 
The lower left panel shows the kinetic temperature we calculate if scattered light is neglected and the lower right panel shows the temperature originally 
computed by KP09. 



3 RESULTS 



3.1 2D simulation (KP09-II) 



Figure 4 shows the mean ionization state of iron and the kinetic 
temperature we compute for the snapshot of the KP09-II model. 



Throughout most of the simulation volume, the gas is very highly 
ionized (the Fe ionization balance is dominated by Fe XXV, XXVI 
and XXVII) and the temperature is very close to the Compton tem- 
perature of our source radiation field (Tc ~ 5 x lO'' K). Only in 
the dense outflowing funnel of gas is the ionization (and tempera- 
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ture) reduced. Even in that region, however, the dominant Fe ion is 
no less ionized than Fe XVII and the temperature remains around 
10^ K, the regime in which the gas is expected to be thermally 
unstable and fragment as found by KP09. Overall, the degree of 
variation in ionization and temperature is significantly less than in 
the disk wind simulation considered by Sim et al. (2010b), a conse- 
quence of the bright primary X-ray source and the absence of any 
very Compton-thick lines of sight through the KP09 simulations. 

Compared to KP09 (see Figure 4), our temperatures are lower 
in the regions of inf ailing gas for two reasons. First, adiabatic 
heating is included in the KP09 simulations but neglected here. 
Moreover, our different assumptions about the spectrum of the pri- 
mary radiation source mean that the Compton temperature in KP09 
(Tc - 2 X 10'' K) is higher than ours (T^ ~ 5 x 10® K). Thus, our 
calculations largely agree that most of the volume of infalling gas 
is heated to around the Compton temperature and the differences 
are mostly due to the assumed source spectral energy distribution. 
This difference does not significantly impact our conclusions since 
even at our Compton temperature, the infalling gas is too hot and 
ionized to imprint strong spectral features (see below). 

In the dense clouds embedded in the outflowing gas we find 
higher temperatures than KP09. Additional radiative heating oc- 
curs in our simulations because of scattered light, which was not 
treated by KP09. To quantify the effect of scattered light, we car- 
ried out an additional calculation in which scattered light was not 
included - the temperature distribution obtained is shown in Fig- 
ure 4. Although less crucial here than in our disk wind simulations 
(Sim et al. 2010b), we find that scattered light is still important in 
heating and ionizing gas that is partially shielded from direct irradi- 
ation by the central source. Further hydrodynamical simulations are 
needed to quantify how this influences the densities and velocities 
of the outflowing gas (see Section 5 for further discussion). 

We have computed spectra for the KP09-II model for a set of 
twenty observer orientations. These uniformly sample the range of 
inclinations < < 1, where 6 = cos~^ /x is the angle between 
the observer's line-of-sight and the polar axis of the model. A sub- 
set of these spectra are shown in Figure 5. 

The synthetic spectra computed for the KP09-II model 
roughly divide into two classes, distinguished by the the observer 
orientation relative to the outflowing funnel of gas. As expected, the 
most complex spectra are predicted for observer orientations closer 
to the pole than the outer boundary of the funnel {6 < 35 deg). Fig- 
ure 6 show an example of this type of spectrum in detail. Roughly 
speaking, the spectrum consists of two components: the direct spec- 
trum of photons from the the primary X-ray source and a contri- 
bution of scattered light generated throughout the simulation vol- 
ume (these components are shown separately in the top left panel 
of Figure 6). The relatively dense gas in the polar funnel imprints 
strong photoelectric and line absorption features in the direct spec- 
trum of the primary source. In contrast, the spectrum of scattered 
light shows only relatively weak absorption since it is predomi- 
nantly formed in the very high ionization gas located outside the 
polar funnel. However, scattering and fluorescence in the outflow 
do give rise to a range of weak emission features, particularly in 
the Fe K region. 

For higher inclination angles {6 > 35 deg), the gas in the ob- 
server line-of-sight is too highly ionized for photoelectric absorp- 
tion to affect the spectrum significantly (see Figure 7). For these 
orientation, the strongest spectral features are the emission lines 
produced by fluorescence in the polar funnel region. In particular, 
there is noticeable net emission in the Fe K region (~ 6.5 keV, a 
blend of individual lines from the K-shell and high L-shell ions). 



yii = cos 6 
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Figure 5. Computed spectra for the KP09-II model. Each line represents a 
different observer inclination angle {d = cos~^ /x) to the polar axis. All 
spectra are normalized to the primary X-ray power-law spectrum and are 
shown at a spectral resolution c HE/E = 5000 km s~^. For clarity, the 
spectra are offset from each other by additive constants (the = 1 spectrum 
is shown with no offset; each of the others is offset by -l-O.l, relative to the 
spectrum for the next higher ^i-value). 

There are a few very weak absorption features in the spectrum, as- 
sociated with the very highly-ionized gas that dominates the vol- 
ume outside the polar funnel. In particular, the Fe XXVI Ka line 
appears as a clear inverse (but weak) P Cygni profile (see lower 
right panel of Figure 6), arising from the net inflow in most of the 
computational grid. Even weaker, but similar, features are also pre- 
dicted in e.g. S XVI Kq and Si XIV Ka (lower left panel of Fig- 
ure 7). 

Since the flow velocities in the KP09 simulations are fairly 
small (at most a few thousand km s~^ but more typically a few 
hundred km s~^), all the Fe K line features in the model spectra 
are relatively narrow such that details of the profiles are not well- 
resolved by current instrumentation (e.g. the Chandra High Energy 
Transmission Grating has a resolution fuU-width-at-half-maximum 
of ~ 1800 km s~^ around the Fe K region). However, the pho- 
toelectric absorption due to the polar funnel material (predicted at 
low inclination) and the net Fe Ka emission (strongest at moderate 
to high inclination) should be detectable. The absence of strong ab- 
sorption signatures associated with the infalling gas is an important 
result. It implies that that the observed lack of red-shifted absorp- 
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Figure 6. Computed spectra for the KP09-II model for an observer inclination angle of fi = cos 9 = 0.925. In the top left panel, the black line shows the 
synthetic spectrum. The red and blue lines show the contributions from photons reaching the observer directly from the primary X-ray source and following 
scattering/reprocessing in the flow, respectively. The remaining panels show details of the spectral regions around the O VIII, Si XIV and Fe Ka features. The 
spectrum is normalized to the primary X-ray power-law spectrum and shown at a spectral resolution corresponding to c AE/E = 1000 km s~^. Monte 
Carlo noise is responsible for percentage-level fluctuations in the spectra. 



tion does not exclude the inflow-outflow pattern of the KP09 (and 
similar) simulations. 



3.2 3D model (KP09-IV) 

The polar symmetry imposed in the 2D model inhibits the devel- 
opment of substructure in the polar funnel. As noted in Section 2, 
in the 3D simulation (KP09-1V) the material in the funnel region 
fragments into clumps and clouds, leading to significant variations 
in the column density through the flow for different azimuthal ori- 
entations (see Figure 3). Here our goal is to quantify how this addi- 
tional variation affects the spectral features. 

To this end, we computed a total of 100 synthetic spectra for 
model KP09-IV. These explore five values of the azimuthal angle 
(0 — 30, 60, 90, 120 and 150 degrees) and the same set of twenty 
^-values used for the KP09-1I spectra. These KP09-IV spectra are 
generally similar to those calculated from KP09-1I and can again 
be subdivided into two broad categories: spectra with strong photo- 
electric absorption (for observer orientations relatively close to the 
pole) and those in which the only significant spectral features are 
narrow fluorescence emission lines. 

As expected (see Figure 3), the 3D flow structure significantly 
increases the diversity of absorption features that can appear in the 



spectrum for near polar inclinations. This is illustrated in Figure 8, 
which compares the synthetic spectra from the KP09-IV simulation 
to the KP09-II model at /i = 0.925. However, for observer orienta- 
tions away from the polar axis (i.e. cases in which there is no strong 
photoelectric absorption), the differences between the KP09-I1 and 
KP09-IV spectra are minor (right panel of Figure 8), indicating that 
the extra structure present in the 3D simulation does not strongly 
influence the fluorescent emission from the polar region. This is ex- 
pected since the emission is formed as an integral over the complete 
flow and neither the total mass outflow rate not the mean volume 
occupied by the flow are very different between the KP09-1I and 
KP09-1V simulations. 



4 DISCUSSION 

Our calculations shows that the KP09 models predict both absorp- 
tion and emission features that could be detectable in X-ray spectra. 
Below, we comment on the comparison of these features with ob- 
servations of AGN. 
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4.1 Absorption features 

Most of the spectral features are associated with the outflowing po- 
lar funnel of gas in the KP09 simulations. For observer orientations 
where the funnel is close to the line-of-sight (9 < 35 deg.), it causes 
significant photoelectric absorption in the 1-10 keV band and also 
imprints a forest of narrow, low- velocity absorption lines associated 
with moderate-to-high ionization gas. 

The absorption components predicted in our simulations are 
not expected to be strongly variable on timescales typically probed 
by X-ray observatories - the simulations cover a region extending 
from ~ 0.01 to 7 pc and have gas velocities of at most a few thou- 
sand km s~^. Therefore the timescale for significant variations is 
expected to be greater than the decade-long timescale of modem X- 
ray missions. Some degree of variation in the degree of ionization 
may be possible (driven by changes in the primary X-ray source 
properties) but this is expected to be small and washed out by the 
long recombination times associated with low-density material at 
large distances (see e.g. Kaastra et al. 2012). 

Therefore, the absorption predicted in our simulations is un- 
likely to be physically related to observed absorption components 
that manifest rapid variability. However, it may be related to per- 
sistently observed absorption systems that do not show strong vari- 
ability and are consequently thought to be relatively distant from 
the AGN. Such absorption components have been identified in a 
number of sources. For example, Kaastra et al. (2012) used the lack 
of short-term variability to constrain the location of several ionized 



absorption components in the X-ray spectra of Mrk 509 to be > 5pc 
from the central black hole. Similarly, Rrongold et al. (2010) con- 
strained the gas responsible for a highly-ionized absorption compo- 
nent in NGC 5548 to lie at > O.OSpc while Behar et al. (2003) and 
Netzer et al. (2003) favour distances of > 0.18 — lOpc for the main 
warm absorbers in NGC 3783 (although we note that Reeves et al. 
2004 argue that the most highly-ionized absorption in NGC 3783 
may be located at smaller radii, ^0.1 pc). These distance limits 
are roughly consistent with the location of the gas in the KP09 sim- 
ulations. Also, the observationally inferred outflow velocities are 
comparable to those in the theoretical models ( < 1000 km s^^). 

However, the total column density associated with the ionized 
absorption by the polar funnels in our simulations and the degree of 
ionization are typically higher than that inferred from the observa- 
tions. For Mrk 509, Kaastra et al. (2012) measure hydrogen column 
densities for their most highly ionized absorption components (D 
and E) of only ~ 6 x 10^" cm"^ and 60 x 10^'' cm~^ while val- 
ues up to ~ lO^'^ cm^^ have been suggested for both NGC 3783 
and NGC 5548 (Behar et al. 2003; Netzer et al. 2003; Krongold 
et al. 2010). In contrast, the KP09 column densities are typically 
~ 10'^'^ cm~^ in the polar region, making the corresponding com- 
ponent of ionized absorption too strong. In addition, observations 
show clear evidence of significant absorption by lower ionization 
gas than is present in the simulations (e.g. the unresolved transition 
arrays of M-shell Fe ions; see Behar et al. 2001). These discrepan- 
cies clearly point to quantitative shortcomings of the properties of 
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the outflow gas in the simulations and motivate further exploration 
of the model parameter space (see Section 5 for further discussion). 

Our calculations predict no strong absorption features associ- 
ated with the infalling gas, despite the significant column densities 
associated with this component of the flow (see Figure 2). Thus, the 
observed absence of strong red-shifted absorption features in X-ray 
spectra does not exclude models for outflows driven from inflows. 
Apart from the weak Fe XXVI inverse P Cygni line (Section 3.1), 
the main effect of the infalling gas is to generating a smooth con- 
tinuum of scattered light (see Section 4.3). 



4.2 Emission features 

Scattering and fluorescence in the outflowing gas gives rise to emis- 
sion lines, including significant Fe Ka. These emission features ap- 
pear for all observer orientations and are the only strong feature for 
inclination angles larger than the opening angle of the polar fun- 
nel. The equivalent width of the Fe Ka emission is ^ 100 eV for 
most observer orientations. This is comparable to the strength of 
the narrow Ka line cores observed in Seyfert galaxies (Yaqoob & 
Padmanabhan 2004; Shu et al. 2010), suggesting that the amount 
of fluorescing gas in the KP09 models may be roughly appropriate. 
However, the position and structure of the synthetic line features 
are not consistent with typically observed narrow lines. In partic- 
ular, Yaqoob & Padmanabhan (2004) and Shu et al. (2010) find 
lines with energy close to 6.4 keV and full-width-at-half-maximum 
(FWHM) of ~ 2000 km s"\ Although the velocities in our mod- 
els are roughly consistent with this line width ( < 1000 km s~^), 
our synthetic spectra contain multiple line components due to the 
range of ionization states present in the outflowing gas (see Fig- 
ures 6 and 7). These different line components could be marginally 
resolved in high-quality data and cause the Fe Ka emission com- 
plex to extend over a much wider spectral region (equivalent to 



~ 10,000 km s~^) than is typically observed. Furthermore, the 
flux is concentrated around 6.5 - 6.6 keV rather than 6.4 keV. Thus, 
compared to the commonly observed properties of narrow Fe K 
emission, it is clear that our calculations produce too much emis- 
sion from significantly ionized material (Fe XIX - XXIV). 

Given that our calculations predict that the Fe Ka emission 
should appear for the majority of observer orientations, the dis- 
crepancies in spectral structure in the region around 6.4 - 6.7 keV 
suggest that the conditions predicted in our current simulations are 
not typical of Type I AGN. However, we note that the discrepan- 
cies could be largely removed if the typical ionization state of the 
material in the outflowing gas were reduced. To qualitatively illus- 
trate the sensitivity of the Fe Ka emission to the flow conditions, 
we have also computed synthetic spectra for the 2D KP09-II model 
in which we fixed the temperature to that obtained by KP09 (as 
discussed in Section 3.1, they generally found lower temperatures 
in the dense regions of the outflow, which result in reduced ioniza- 
tion compared to our full radiative transfer calculation). The syn- 
thetic spectrum from this calculation is compared to our standard 
case in Figure 9. This shows that the Ka emission is affected by 
details of the conditions in the outflow - adopting the KP09 tem- 
peratures changes the ionization state sufficiently that more of the 
flux emerges around 6.4 keV. However, this change alone does not 
qualitatively alter our findings - a more significant reduction of the 
ionization state would be required to shift most of the Ka emission 
to 6.4 keV. This might be achievable via stronger clumping in the 
outflow, as will be discussed in Section 5. 

We note that our models also predict an emission compo- 
nent from Fe XXVI around 6.9 keV, which is well-separated from 
the main Fe Ka emission in our spectra (6.4 - 6.7 keV). Narrow 
Fe XXVI emission components have been reported in several ob- 
jects (e.g. NGC 7314 [Yaqoob et al. 2003]; NGC 7213 [Bianchi 
et al. 2008; NGC 3516 [Turner et al. 2008]) and based on a sample 
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of Suzaku spectra, Patrick et al. (20 11) suggested ttiat such a feature 
is relatively common (although we note that this result is somewhat 
dependent on the modelling, see Tatum et al. 2012). 

Particularly for inclination angles beyond the polar funnel 
(6 > 35 deg.), our synthetic spectra also contain emission lines as- 
sociated with H-like ions of lighter elements, most prominently 
O VIII Lyman a (but also S XVI, Si XIV and Mg XII Lyman a). 
Such H-like emission lines have been reported in e.g. the spec- 
tra of the Seyfert 2 galaxy NGC 1068 (Kinkhabwala et al. 2002; 
Ogle et al. 2003). This is qualitatively consistent with our simula- 
tions, which predict that emission lines are the only strong spectral 
features imprinted for relatively high inclinations. However, there 
are important differences. In particular, the observations also in- 
clude many lines of lower ionization material (e.g. Kinkhabwala 
et al. 2002 and Ogle et al. 2003 report strong emission in both 
O VII and O VIII while the model spectra show virtually no O VII 
emission lines). There is also a kinematic mismatch - e.g. the syn- 
thetic O VIII line profile is not resolved in our Monte Carlo spec- 
tra (FWHM < 500 km s^^) while the observed profile is broader 
(FWHM ~ 1500 km s~^), suggesting that higher fluid velocities 
are required. 

To summarise, comparison of the emission line properties 
leads to similar conclusions to those drawn from the absorption 
features: the characteristic properties of the calculated spectra are 
encouraging but there are clear quantitative shortcomings. In par- 
ticular, the ionization state is generally too high in the simulations. 



4.3 Scattered continuum 

In all our synthetic spectra, a significant fraction of the predicted 
continuum flux (> 10 per cent; see the scattered light component 
in Figures 6 and 7) is produced by Compton scattering in the highly 
ionized gas outside the polar funnel. Consequently, even for ori- 
entations for which the absorbing column through the funnel is 
very high, the observer will see a component of weakly attenuated 
continuum flux (produced by Compton scattering into the line-of- 
sight). As a result, absorption features are not expected to reach 
zero flux, despite high optical depths. The spectrum can therefore 
resemble that of a partial-covering model (in which the spectrum is 



described by a linear combination of components, typically one un- 
absorbed continuum component plus a significantly absorbed com- 
ponent). An example of scattered continuum re-filling absorption 
lines, albeit on a different scale (150 pc), has been seen in the case 
of Mrk 509 (Kriss et al. 2011). 



5 SUMMARY AND FUTURE WORK 

We have calculated synthetic X-ray spectra for the KP09 two- and 
three-dimensional simulations of outflows driven from infalling gas 
close to an AGN. We showed that the polar funnel of outflowing gas 
predicted by the KP09 simulations has observable consequences for 
all observer orientations. It causes strong, ionized absorption fea- 
tures for near face-on inclinations and ionized fluorescent emission 
for all lines of sight. The shapes of the spectral features predicted 
for our two- and three-dimensional simulations are generally simi- 
lar. However, in the three-dimensional case, we do find significant 
variation (a factor of ~ 3) in the depth of absorption with azimuthal 
angle (see Figure 8), which is a consequence of the clouds/clumps 
in the flow. 

Conversely, we found that the infalling gas, which occupies 
most of the volume outside the polar funnel in the simulations, does 
not imprint strong spectroscopic absorption features because of the 
very high ionization state. Only a weak Fe XXVI Ka feature is pre- 
dicted and we note that our simulations may even over-estimate the 
strength of this feature (non-radiative heating likely makes the tem- 
perature of the infalling gas even higher; see Section 3.1). Thus we 
conclude that an observed lack of inflow features does not rule out 
outflow-from-inflow models similar to those presented by KP09. 

When compared to the properties of observed absorption and 
emission features associated with gas that is relatively distant 
from the black hole, the particular simulations presented here face 
several important challenges. Although the outflow velocities are 
roughly consistent with observations, the column density and typi- 
cally ionization state of the outflowing gas are too high in the mod- 
els. 

The ionization state in our simulations could be lowered by 
simply reducing the source X-ray luminosity. However, the reduc- 
tion required (a factor of ~ 10) would be quite inconsistent with 
the properties of the radiation source assumed in the KP09 hydro- 
dynamical simulations. Moreover, reduced luminosity leads to a 
global reduction in the ionization of all the gas - consequently, the 
infalling gas (occupying ^ 35 deg) would cease to be fully ion- 
ized and start to imprint spectral features. The fact that the absorp- 
tion predicted for polar orientations {9 < 35 deg) is also too strong 
in the simulations makes this difficulty harder to overcome. If the 
flow pattern remains the same, the gas density would need to be 
lower to give a smaller column density but a reduced density would 
make the over-ionization problem more severe. 

More promising than any global change of the density or X- 
ray flux might be smaller scale changes in the structure of the out- 
flowing gas in the polar funnel. For a similar total mass outflow rate, 
higher typical densities in the outflow might be achieved if it were 
more strongly clumped. This would not only lead to lower ioniza- 
tion but could also reduce the column density since smaller clumps 
will subtend smaller solid angles as seen by the central source. 

It is important to note that the KP09 simulations, and the as- 
sociated radiative transfer calculations presented here, are only ini- 
tial, exploratory results. A major conclusion to be drawn is that 
thermal instabilities are likely to cause the flow to fragment, which 
has observable consequences - that the thermal instability may be 
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at work in AGN outflows is already supported by spectroscopic ev- 
idence (Holczer et al. 2007). More detailed work is now needed 
to explore the parameter space of the simulations and investigate 
physical and numerical aspects of the fragmentation process. For 
example, here we have considered only one case for the angular 
momentum of the infalling gas (corresponding to a circularization 
radius of rc = 1800r-g). Changing the angular momentum will al- 
ter the relative densities of polar and equatorial parts of the flow. 
Moreover, although the KP09 simulations do show that the out- 
flowing gas breaks into clouds, the scale of fragmentation is lim- 
ited by the resolution of the simulations and depends on factors 
such as the accuracy of the opacities used in the hydrodynamical 
simulations. Further investigation of fragmentation in simulations 
for gas infalling towards an AGN have been made by Barai et al. 
(201 1, 2012). They carried out simulations extending to larger spa- 
tial scales than KP09 but did not include the radiative line-force. 
Their results confirm that thermal instabilities can cause clumps to 
form and, depending on simulation parameters, that this can also 
happen in the infalling gas. How such fragmentation of the inflow 
might affect the dynamics of a radiatively driven outflow and the 
observables now needs to be investigated via simulations similar to 
KP09 with extended computational domains. In addition, our ra- 
diative transfer simulations have demonstrated that scattered light 
can be significant in heating and ionizing the outflowing material 
(see also Sim et al. 2010b), an effect which is not currently incor- 
porated in the hydrodynamical simulations. Simply adopting the 
temperature distribution derived by KP09 in our radiative trans- 
fer simulations does not qualitatively change our results (although 
our synthetic spectra are quantitatively affected). Nevertheless, hy- 
drodynamical simulations that consider scattered light are needed 
since the additional heating may affect the onset of thermal insta- 
bilities and ultimately impact the density and velocity distribution 
of the outflowing gas. 

We conclude that outflows driven from inflows are promising 
as a possible site for observable X-ray absorption and fluorescent 
emission around AGN. However, further study is required: in par- 
ticular, the next step should be detailed, high-resolution studies of 
the clumping and structure of outflows and how these depend on 
key simulation parameters. 
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